Search for: All records

Abstract Inland waters receive large quantities of dissolved organic carbon (DOC) from soils and act as conduits for the lateral transport of this terrestrially derived carbon, ultimately storing, mineralizing, or delivering it to oceans. The lateral DOC flux plays a crucial role in the global carbon cycle, and numerous models have been developed to estimate the DOC export from different landscapes. We reviewed 34 published models and compared their characteristics to identify challenges in model applications and opportunities for future model development. We classified these models into three types: indicator-driven, hydrology-forced, and process-based DOC export simulation models. They differ mainly in their environmental inputs, simulation approaches for soil DOC production, leaching from soils to inland waters, and transit through inland waters. It is essential to consider landscape characteristics, climate conditions, available data, and research questions when selecting the most appropriate model. Given the substantial assumptions associated with these models, sufficient measurements are required to benchmark estimates. Accurate accounting of terrestrially derived DOC export to oceans requires incorporating the DOC produced in aquatic ecosystems and deposited with rainwater; otherwise, global export estimates may be overestimated by 40.7%. Additionally, improving the representation of mineralization and burial processes in inland waters allows for more accurate accounting of carbon sequestration through land ecosystems. When all the inland water processes are ignored or assuming DOC leaching is equivalent to DOC export, the loss of soil carbon through this lateral flux could be underestimated by 43.9%. 

The positive Arctic–methane (CH₄) feedback forms when more CH₄is released from the Arctic tundra to warm the climate, further stimulating the Arctic to emit CH₄. This study utilized the CLM-Microbe model to project CH₄emissions across five distinct Arctic tundra ecosystems on the Alaska North Slope, considering three Shared Socioeconomic Pathway (SSP) scenarios using climate data from three climate models from 2016 to 2100. Employing a hyper-resolution of 5 m × 5 m within 40,000 m²domains accounted for the Arctic tundra’s high spatial heterogeneity; three sites were near Utqiaġvik (US-Beo, US-Bes, and US-Brw), with one each in Atqasuk (US-Atq) and Ivotuk (US-Ivo). Simulated CH₄emissions substantially increased by a factor of 5.3 to 7.5 under the SSP5–8.5 scenario compared to the SSP1–2.6 and SSP2–4.5 scenarios. The projected CH₄emissions exhibited a stronger response to rising temperature under the SSP5–8.5 scenario than under the SSP1–2.6 and SSP2–4.5 scenarios, primarily due to strong temperature dependence and the enhanced precipitation-induced expansion of anoxic conditions that promoted methanogenesis. The CH₄transport via ebullition and plant-mediated transport is projected to increase under all three SSP scenarios, and ebullition dominated CH₄transport by 2100 across five sites. Projected CH₄emissions varied in temperature sensitivity, with a Q₁₀range of 2.7 to 60.9 under SSP1–2.6, 3.8 to 17.6 under SSP2–4.5, and 5.7 to 17.2 under SSP5–8.5. Compared with the other three sites, US-Atq and US-Ivo were estimated to have greater increases in CH₄emissions due to warmer temperatures and higher precipitation. The fact that warmer sites and warmer climate scenarios had higher CH₄emissions suggests an intensified positive Arctic–CH₄feedback in the 21st century. Microbial physiology and substrate availability dominated the enhanced CH₄production. The simulated intensified positive feedback underscores the urgent need for a more mechanistic understanding of CH₄dynamics and the development of strategies to mitigate CH₄across the Arctic. 

Abstract Large across‐model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across‐model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO 2 levels and five temperature levels. We developed a matrix‐based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter‐model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter‐model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO 2 (eCO 2 ) concentrations in the measurements, simulated NEP responded positively to eCO 2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across‐model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables. 

Spatial heterogeneity in methane (CH 4 ) flux requires a reliable upscaling approach to reach accurate regional CH 4 budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH 4 flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH 4 flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH 4 flux in early growing seasons. The simulated monthly CH 4 flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m −2  s −1 and 0.782 using the DF algorithm, but 0.007 μmol m −2  s −1 and 0.758 using HF and 0.007 μmol m −2  s −1 and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH 4 flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH 4 flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH 4 flux for all five tower domains despite relatively weak differences in simulated CH 4 flux among three footprint algorithms. The CLM-Microbe model can simulate CH 4 flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH 4 flux in terrestrial ecosystems. 

Urbanization has caused environmental changes, such as urban heat islands (UHIs), that affect terrestrial ecosystems. However, how and to what extent urbanization affects plant phenology remains relatively unexplored. Here, we investigated the changes in the satellite-derived start of season (SOS) and the covariation between SOS and temperature ( R T ) in 85 large cities across the conterminous United States for the period 2001–2014. We found that 1) the SOS came significantly earlier (6.1 ± 6.3 d) in 74 cities and R T was significantly weaker (0.03 ± 0.07) in 43 cities when compared with their surrounding rural areas ( P < 0.05); 2) the decreased magnitude in R T mainly occurred in cities in relatively cold regions with an annual mean temperature <17.3 °C (e.g., Minnesota, Michigan, and Pennsylvania); and 3) the magnitude of urban−rural difference in both SOS and R T was primarily correlated with the intensity of UHI. Simulations of two phenology models further suggested that more and faster heat accumulation contributed to the earlier SOS, while a decrease in required chilling led to a decline in R T magnitude in urban areas. These findings provide observational evidence of a reduced covariation between temperature and SOS in major US cities, implying the response of spring phenology to warming conditions in nonurban environments may decline in the warming future. 

Abstract Satellite‐derived sun‐induced chlorophyll fluorescence (SIF) has been increasingly used for estimating gross primary production (GPP). However, the relationship between SIF and GPP has not been well defined, impeding the translation of satellite observed SIF to GPP. Previous studies have generally assumed a linear relationship between SIF and GPP at daily and longer time scales, but support for this assumption is lacking. Here, we used the GPP/SIF ratio to investigate seasonal variations in the relationship between SIF and GPP over the Northern Hemisphere (NH). Based on multiple SIF products and MODIS and FLUXCOM GPP data, we found strong seasonal hump‐shaped patterns for the GPP/SIF ratio over northern latitudes, with higher values in the summer than in the spring or autumn. This hump‐shaped GPP/SIF seasonal variation was confirmed by examining different SIF products and was evident for most vegetation types except evergreen broadleaf forests. The seasonal amplitude of the GPP/SIF ratio decreased from the boreal/arctic region to drylands and the tropics. For most of the NH, the lowest GPP/SIF values occurred in October or September, while the maximum GPP/SIF values were evident in June and July. The most pronounced seasonal amplitude of GPP/SIF occurred in intermediate temperature and precipitation ranges. GPP/SIF was positively related to temperature in the early and late parts of the growing season, but not during the peak growing months. These shifting relationships between temperature and GPP/SIF across different months appeared to play a key role in the seasonal dynamics of GPP/SIF. Several mechanisms may explain the patterns we observed, and future research encompassing a broad range of climate and vegetation settings is needed to improve our understanding of the spatial and temporal relationships between SIF and GPP. Nonetheless, the strong seasonal variation in GPP/SIF we identified highlights the importance of incorporating this behavior into SIF‐based GPP estimations. 

Abstract Transitioning across biological scales is a central challenge in land surface models. Processes that operate at the scale of individual leaves must be scaled to canopies, and this is done using dedicated submodels. Here, we focus on a submodel that prescribes how light and nitrogen are distributed through plant canopies. We found a mathematical inconsistency in a submodel implemented in the Community and Energy Land Models (CLM and ELM), which incorporates twigs, branches, stems, and dead leaves in nitrogen scaling from leaf to canopy. The inconsistency leads to unrealistic (physically impossible) values of the nitrogen scaling coefficient. The mathematical inconsistency is a general mistake, that is, would occur in any model adopting this particular submodel. We resolve the inconsistency by allowing distinct profiles of stems and branches versus living leaves. We implemented the updated scheme in the ELM and find that the correction reduces global mean gross primary production (GPP) by 3.9 Pg C (3%). Further, when stems and branches are removed from the canopy in the updated model (akin to models that ignore shading from stems), global GPP increases by 4.1 Pg C (3.2%), because of reduced shading. Hence, models that entirely ignore stem shading also introduce errors in the global spatial distribution of GPP estimates, with a strong signal in the tropics, increasing GPP there by over 200 g C m⁻² yr⁻¹. Appropriately incorporating stems and other nonphotosynthesizing material into the light and nitrogen scaling routines of global land models, will improve their biological realism and accuracy.